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Thermodynamics of the six-vertex model in an L-shaped 

domain 

Filippo Colomo and Andrei G. Pronko 


Abstract. We consider the six-vertex model in an L-shaped domain of the 
square lattice, with domain wall boundary conditions. For free-fermion ver¬ 
tex weights the partition function can be expressed in terms of some Hankel 
determinant, or equivalently as a Coulomb gas with discrete measure and a 
non-polynomial potential with two hard walls. We use Coulomb gas methods 
to study the partition function in the thermodynamic limit. We obtain the 
free energy of the six-vertex model as a function of the parameters describing 
the geometry of the scaled L-shaped domain. Under variations of these pa¬ 
rameters the system undergoes a third-order phase transition. The result can 
also be considered in the context of dimer models, for the perfect matchings 
of the Aztec diamond graph with a cut-off corner. 


1. Introduction 

It is commonly known that macroscopic quantities (such as the free energy) of 
dimer coverings and random tilings on regular lattices may depend on boundary 
conditions m- This feature is related to phase separation and the emergence of a 
limit shape [3H7] , see [5] for a review. The same phenomena can be observed in the 
six-vertex model [SHII]. In relation with these phenomena, an interesting question 
concerns the stability of the observed bulk properties against various deformations 
of the shape of the considered finite region, while preserving the particular boundary 
conditions (e.g., the staircase shape for the boundary of the Aztec diamond in the 
case of domino tilings). 

In this paper, we address the problem by studying the thermodynamics of the 
six-vertex model on a particular domain of the square lattice, that we call an L- 
shaped domain. The model is closely related to the domino tilings of the Aztec 
diamond with a cuf-off corner [12]. We evaluate the exact analytic expression for 
the free energy per site of the six-vertex model at its free fermion point, as a 
function of the parameters describing the geometry of the scaled L-shaped domain. 
Under variation of these parameters, when the boundary interferes with the phase 
separation curve, the system undergoes a third-order phase transition. 

Our result on the free energy of the six-vertex model in an L-shaped domain 
is based on some discrete Coulomb gas representation. Discrete Coulomb gases 
already appeared in applications to tiling and dimer models HKH], and to the six- 
vertex model [laiisj, and have been studied intensively in the context of discrete 
orthogonal polynomials, see m and references therein. As a side result, we find 
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Figure 1. (a) The L-shaped domain with domain wall boundary 
conditions; (b) The corresponding arrow configuration on the N x 
N lattice. 


that the soft-edge/hard-edge transition may induce (or not) a third-order phase 
transition in the free energy, according to the finite (or diverging) slope of the 
discrete Coulomb gas potential in the vicinity of the hard wall. 

1.1. The model. The states of the six-vertex model are configurations of 
arrows on the edges of a square lattice, satisfying the ice rule: at each vertex the 
numbers of incoming and outgoing arrows are equal. The ice rule selects ( 2 ) = 6 
possible vertex configurations, hence the name of the model. 

The L-shaped domain can be defined as a square domain with a rectangular 
portion removed from one of the corners, see Fig.[TJ Specifically, the square domain 
is the finite square lattice obtained from the intersection of N horizontal and N 
vertical lines, the so-called N x N lattice. The L-shaped domain is obtained by 
removing from the top left corner a rectangular portion of the lattice, of size s x 
(s -I- q), s,q G No, s + q < N. The L-shaped domain is shown in Fig. [T^, where 
r = N — s — q. 

We consider a specific instance of fixed boundary conditions, the domain wall 
boundary conditions. In the standard arrows picture of the six-vertex model (see, 
e.g., [18]) this conditions means that all horizontal (respectively, vertical) arrows 
on external edges are outgoing (incoming). For s = 0, one has the usual N x N 
lattice with domain wall boundary conditions, introduced in |19j . 

The Boltzmann weights of the six-vertex model, Wi, i = 1,... ,6, enumerated 
in the standard order, see Fig. O second row, are chosen to satisfy 

W1W2 + W3W4 = W5WQ, ( 1 . 1 ) 

that is called the free-fermion condition. 

We denote the partition function of the six-vertex model in an L-shaped domain 
with domain wall boundary conditions by Zr^s,q] for the special case s = 0, we use 
the standard notation Zjq. For generic Boltzmann weights, the partition function 
Zjq has been evaluated in determinantal form in |20l[2T] . For weights satisfying the 
condition (EB , Z]\[ has an extremely concise form: 

iV(iV-l) N{N + 1) 

Zn =w^ ^ Wq ^ . 
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Figure 2. Small patches of dominoes (first row) and vertex con¬ 
figurations of the six-vertex model (second row). 




(a) (b) 

Figure 3. (a) Aztec diamond with a frozen region of NE-SW 
dominoes; (b) Aztec diamond with a cut-off corner and the corre¬ 
sponding modified Aztec diamond graph. 


For a proof, see [5] in the context of enumerative combinatorics, or |22] in the 
context of integrable models. 

For generic s, the partition function Zr,s,q can be expressed in terms of certain 
non local correlation function of the model on the N x N lattice. This function is 
called the emptiness formation probability, and denoted here as Fr^s,q- It describes 
the probability that all vertices in the top left sx (s+q) rectangle of the NxN lattice 
are of type 2. The corresponding configuration is shown in Fig.jTjD. Comparing the 
two pictures of Fig. [U it is clear that 


Zrcfi - 


F 


( 1 . 2 ) 


Thus the partition function of the six-vertex model on the L-shaped domain is 
essentially given by the emptiness formation probability of the model on the NxN 
lattice. 

The above construction can be translated in the language of dimer models, 
using the well known correspondence between the six-vertex model with domain 
wall boundary conditions and the domino tilings of the Aztec diamond [2] . 

Recall that domino tilings are coverings of a square lattice with rectangles of 
size 1x2 and 2x1. Domino tilings can be formulated in terms of the six-vertex 
model on a square lattice by mapping elementary patches of domino tilings to 
arrow configurations as shown in Fig. [51 The case of the six-vertex model on an 
NxN square lattice with domain wall boundary conditions corresponds to the 
Aztec diamond of order N, see Fig. O 

The plain enumeration of domino tilings of the Aztec diamond corresponds to 
the choice of the Boltzmann weights of the six-vertex model wi = ■ ■ ■ = = 
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1 and wq = 2. In counting of domino tilings, one may also consider weighted 
enumerations, with a bias parameter a € [0,1], describing the asymmetry between 
the two possible orientations of dominoes. Following paper [5], we assign to the 
NE-SW dominoes the weight y^2(l — a), and to the NW-SE the weight To 

establish the correspondence with the six-vertex model, it is useful to parameterize 
the Boltzmann weights as follows: 

Wi = W2 = \/p(l - a), W3 = W4 = y/pa, W5 = 1, We = p. (1.3) 

Then the weighted enumeration of dominoes corresponds to setting p = 2. 

In the same spirit, the six-vertex model in an L-shaped domain can be related 
to the domino tilings of some suitable modification of the Aztec diamond. Using 
the correspondence of vertices of type 2 with a particular patch of domino tilings, 
see Fig. [21 it is clear that emptiness formation probability is equivalent to the 
probability of observing the domino configuration shown in Fig. [3^. Removing the 
frozen region of NE-SW dominoes from the Aztec diamond, one obtain the Aztec 
diamond with a cut-off corner shown in Fig. |3 }d. The same pictures also shows the 
resulting graph for dimer coverings (modified Aztec diamond graph). 


1.2. Discrete Coulomb gas representation. As discussed above, the par¬ 
tition function of the six-vertex model on an L-shaped domain is essentially given 
by the emptiness formation probability. For generic values of Boltzmann weights, 
the emptiness formation probability has been computed in the form of a multiple 
integral in |23j . Under the free-fermion condition this representation reduces to an 
Hankel determinant [24] . 


Proposition 1.1. The emptiness formation probability of the six-vertex model 
on the N X N lattice with domain wall boundary condition, with the Boltzmann 
weights (lESl admits the following representation: 


p = 


{q'r 


(1-a) 


Ks+'j) 


nrJo(<i+ky-k'- 


ySis-l)/2 


det 

l<j,k<s 



(1.4) 


For a proof, see [231124] . 


Remark 1.2. An alternative representation, equivalent to (d, reads 

r, (1 - r 


where 


^r,s,q — 


in 


q'- 


“ q,sU-1)/2 

E-E n irn,-m,rfl 


s’ fhf-f-oV 

mi=0 m^=01<j<k<s 


i=i 


q ruj 
q 


(1.5) 


( 1 . 6 ) 


The last representation is that of a discrete Coulomb gas confined within a 
finite interval. We recognize in the discrete weight 

Wq{m) = m e No, 

that of Meixner polynomials, but we emphasize the condition m < r for the charge 
coordinates in the Coulomb gas representation (ITBl) . This condition plays a major 
role in what follows. We shall refer to this condition as an hard wall for the charges. 
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The discrete Coulomb gas formula (HH) appeared previously in the context of 
some random growth model studied by Johansson m- The connections of (ira 
with the circular unitary ensemble, with a multivariate generalization of hyper¬ 
geometric function, and with the r-function of the sixth Painleve equation were 
enlightened by Forrester and Witte [25] . 

1.3. Thermodynamics: the main result. We are interested in the be¬ 
haviour of the model in the thermodynamic limit, i.e., in the limit of the large 
lattice sizes, r,s,q —>■ oo. This is described by the free energy per site, defined as 
follows. Let s + q = [xA^], s = [yN], with x,y G [0,1] fixed, x > y. We set 

/(a;,y) := - lim —^— log (1-7) 

N^oo — s{s + q) 

Equivalently, we may study the emptiness formation probability Fr^s.q- We set 

(p{x,y) := - \im log FN-s-q,s,q- (1-8) 

N^oo 

Remark 1.3. The above limits exist; in particular, ip{x,y) is non-negative all 
over its domain of definition, x,y G [0,1], x> y. 

This is part of a stronger statement proven in m, see Theorem 1.1 therein. 

From relation m it immediately follows that the free energy per site of the six- 
vertex model on the L-shaped domain is completely determined by the knowledge 
of the function ip{x,y). Indeed, one has 

(1 - xy)f{x, y) = /(0,0) -b xy log W 2 + y}{x, y), /(0,0) = - log 

Similarly, the function (p{x, y) determines completely the free energy of the domino 
tilings (up to trivial modifications, which can be easily inferred from (|1.2ll L 

The main resnlt of the present paper concerns the explicit form of the function 
T{x,y)- 

Let us consider the unit square, parameterized hy x,y G [0,1]. It follows from 
m that the domain of definition of the function ip{x, y) is the triangular region 
delimited by the three lines y = 0, y = x, and y = 1 — a;. Let us consider the ellipse 
defined by the equation 

(1 — a; — y)^ (x — y)^ 

^^ -b L.;-^ = 1. 1.9 

a 1 — a 

We denote by Vi and Vn the two regions of the triangle that are inside, respectively 
outside, the ellipse, see Fig. IH The arc of the ellipse (jl.9ll lying in the triangle, and 
separating these two regions, is described by the equation: 

y/y = \/{1 — x){l — a) — \fxa, D < y < x <1 — a. (1-10) 

We denote this arc by A. 

Theorem 1.4. The function (p{x,y), {x,y) G [0,1] x [0,1], is 
(p{x,y) = 0, {x,y)GVi, 
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Figure 4. The Arctic ellipse (dotted), the arc A (bold), and the 
domains and Pn. 


and 


h {1-x-y)^ 1-h 1 1 + h 

^{x, y) = xy log - ^13^ - 2 TT^ 

, y + il — x)h , ,, x + [1 — y)h 

+ (1 - x)y log —-— + x{l - y) log ■ 


y + {l- x)r] 

- (1 - x)x log ^"1" ^ - (1 - y)y log 


{x - y)^ 


log 


X + (1 — x)ri 
x + y+(2 — X — y)h 
x + y + {2 — X — y)ri ’ 


a; + (1 - y)v 

y+{l-y)h 
y+{l- y)v 

{x,y) e Vii, 


where 


h = h{x,y) := 

and rj = r]{x, y\ a) is such that 

[ 0 , 1 ], 


xy 


{l-x){l-y) 

(1 + r])^{x + (1 - x)r]) (y + (1 - y)r]) 


(1 - ? 7 ) 2 (?/ + (1 - x)r]) (a: + (1 - y)r)) 


= 1 . 


( 1 . 11 ) 


( 1 . 12 ) 


(1.13) 


The fact that ip{x,y) vanishes in the domain was already proven in |13j . A 
few comments are in order. 


Remark 1.5. The function ip{x,y) is a continuous function of its variables. In 
particular it vanishes on the curve A: 

This property can be verified as follows. Let {x,y) S A, then a = Oc, where 

Oc = ac(x, y) := (^(1 - a;)(l - y) - ■ (1.15) 
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Setting a = ac in (11-13I) . one can directly check that the root of equation (11.131) is 
r] = h, with h given by (11.121) . In other words, 

v{x, y; a) I = y(x, y, ac) = h{x, y), (1.16) 

and (11.141) follows from (11.111) . 

The function y) = <f{x, y; a) may also be viewed as a function of a at fixed 
X and y. 


Remark 1.6. The function (p{x,y) satisfies 

(p{x,y;ac) = 0. 

This property follows from the second equality in (I1.16p . In fact, stronger 
properties hold. 

Proposition 1.7. Given {x,y), for values of a in the vicinity of ac = ac{x,y), 
we have 

(p{x,y) = 

where C = C{x,y) > 0. 


0 a < ac 

C{a - ac)^ + O ({a - ac)^) a > ac, 


Proposition 1.8. Given a point {x,y) € Vu, at a distance t from the arc of 
ellipse A, one has 

<p{x, y) (xt^, t —>■ O’*". 

where the proportionality constant is positive. 


These properties follows from the explicit form of function ip{x, y) given in 
Theorem ll.41 


1.4. Discussion. The ellipse p.9l) is easily recognized as the Arctic ellipse of 
the domino tilings of the original, non-deformed Aztec diamond [3], or, equivalently, 
of the six-vertex model with domain wall boundary condition, at its free-fermion 
point, on the original square domain. 

The function ip{x,y) vanishes at a = ac, together with its first and second 
derivatives with respect to a, and has non-vanishing (and finite) third-order de¬ 
rivative. This property obviously holds also for the free energy f{x,y), see p.7D . 
Recalling that the parameter a is directly related to the Boltzmann weights, we 
thus observe the occurrence of a third-order phase transition at a = ac- 

An alternative point of view is that of considering a fixed, while varying the 
shape of the L-shaped domain, and thus the coordinates x and y. The function 
ip{x,y), studied as a function of the parameter t, again undergoes a third-order 
phase transition. 

This result leads to an alternative interpretation of the Arctic ellipse, as a 
critical curve in the space of the parameters of the model, a point of view already 
advocated in [TOj. The result extends that reported in [12], where the analysis was 
restricted to the line x = y, and only the critical point {x, y) = () could 
be observed. 

The vanishing of the function (p{x, y) in domain is somewhat obvious, since 
it is known that in (Dj, the emptiness formation probability is equal to 1 up to 
exponentially small corrections, which are related to the so-called upper tail rate 
function [13]. The function (p{x, y) in domain Vu defines instead the so-called lower 
tail rate function, whose explicit form is given by (11.111) . 
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In the derivation of Theorem (d we rely heavily on representation (ira . that 
is a Coulomb gas on a discrete lattice, confined by a non-polynomial potential 
with two hard walls. Thus, a by-product of our results is that the hard wall in the 
discrete Coulomb gas can induce a soft-edge/hard-edge transition with a third-order 
discontinuity in the corresponding free energy. This is a well-known phenomenon 
in the case of a continuous Coulomb gas, see [26] and references therein. However, 
in the discrete case, phase transitions are rather of the Douglas-Kazakov type m, 
while the soft-edge/hard-edge transition is usually transparent, in that it does not 
come with any discontinuity in the free energy (see m for numerous examples). 

The peculiar behaviour observed here can be ascribed to the form of the po¬ 
tential in the vicinity of the hard walls. In previously studied cases, in the vicinity 
of a hard wall situated, say, at xq, the derivative of the potential has a divergent 
behaviour V{x) ~ log \x — a;o|i implying a smooth matching of the potential with 
the hard wall. In the present model, instead, such smooth matching occurs only 
in correspondence of the left hard wall. As a result, among the two possible soft- 
edge/hard-edge transitions in the present model only one is transparent, that is 
related to the left hard wall. The other one, related to the right hard wall, comes 
indeed with a third-order phase transition in the corresponding free energy. 

Acknowledgments. We thank A. Kuijlaars and P. Zinn-Justin for useful 
discussions. This work is partially supported by the IRSES grant of EC-FP7 
Marie Curie Action “Quantum Integrability, Conformal Field Theory and Topo¬ 
logical Quantum Computation”. A.G.P. acknowledges partial support from the 
Russian Foundation for Basic Research, grant 13-01-00336, and from INFN, Sezione 
di Firenze. 


2. The equilibrium measure 

We first outline some general aspects of the asymptotic behaviour of the discrete 
Coulomb gas, along the lines developed in see also [T31II1], and next switch 
to the discussion of explicit solutions for the resolvent of the equilibrium measure. 
The end-point equations and the free energy are discussed in the next section. 

2.1. The discrete Coulomb gas. In order to evaluate the limit (HH) and 
prove Theorem 11.41 we study the asymptotic behaviour of the discrete Coulomb 
gas ([T^ in the scaling limit, that is r, s, g —>■ oo, while preserving the aspect ratio 
of the L-shaped domain. Let r = [As], i? > 1 fixed and q = [Qs], Q > 0 fixed. We 
set 

^{R,Q):= lim -^\ogIr,s,q- (2.1) 

s—>-oo 

It was shown in [13], see Theorem 2.2 therein, that the above limit exists. The 
function ^{R,Q) is related to the function (p{x,y) introduced in (11.811 by 

(1 - xfip{x,y) = -xy\og{l - a) -f log Va - (-—, (2.2) 

\ y y J 

see dOI). To obtain relation (12.21) one has to take into account the relations 

R=^, Q = ^^, (2.3) 

y V 

which hold in the limit s —>■ oo. 






V{f^) 


0 



Figure 5. The potential V{fi) with two hard walls, at /i = 0 and 
fj, = R. 

In [12], the focus was on the distribution of the position of the rightmost charge, 
maxi<j<;,mj, that was shown to follow, in the limit s —>■ oo, the Tracy-Widom 
distribution [55] . In m, the explicit form of the upper tail rate function was 
evaluated. Our aim here is to evaluate the explicit form of function $(i?, Q), and 
hence, the lower tail rate function ip{x,y). 

As discussed in m to evaluate the continuous limit of the discrete Coulomb 
gas (d^ we may just as well consider the corresponding continuous Coulomb gas, 
in a continuous potential derived from the original discrete measure. Set 



Using Stirling formula, we have 


V{fi):= \im Vsin) =-fi\oga +ulogfi-{Q + n)log{Q + y) + QlogQ, (2.4) 


uniformly on compact subsets of [0,oo). Note that the confining potential U(/i) is 
complemented by two hard walls, at /r = 0 and /i = i?, see Fig. jS] We emphasize 
the difference in the behaviour of the potential in the vicinity of the two hard walls: 
the derivative of the potential diverges logarithmically at /r = 0, while it stays finite 
at n = R. As a result, the potential joins up the hard wall smoothly on the left, 
but with a corner on the right. 

To proceed, we rely on some known results from weighted potential theory in 
the case of a discrete Coulomb gas. This has been discussed both in the physical [ 57 ] 
and mathematical literature [ 551120 ] ; see for a general mathematical treatment. 
In a nutshell, under mild assumptions on the potential that are clearly satisfied 
in our problem, see |13] . the asymptotic behaviour of the discrete Coulomb gas is 
governed by a constrained equilibrium problem analogous to that for the continuous 
Coulomb gas, with the additional constraint that the particle density is always less 
than or equal to 1, as a consequence of discreteness. 

More specifically, the problem reduces to the evaluation of the equilibrium 
measure, po(/^)d/i, that is the unique minimizer of the functional 


subject to the constraint (normalization condition) 




(2.5) 
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and to the additional constraint 

p{^i) < 1 , ( 2 . 6 ) 

that is imposed by the discreteness of the measure of the Coulomb gas. In the 
above formulas, integration is understood over the interval 

It is known that, due to the presence of the last constraint, the equilibrium 
measure partitions the interval [0, i?] into three sets, , /q, /+, according to po{p) = 
0, 0 < Po(m) < 1) or po{p = 1, respectively. Each set is the union of a finite 
number of intervals, named voids, bands, or saturated regions, respectively. Voids 
and saturated regions are collectively named gaps. 

2.2. The resolvent. The most convenient way to evaluate the equilibrium 
measure is by introducing its resolvent 

W{z)=r^dp, zi[0,R], 

Jo z — p 

that has the following properties: 

(1) W{z) is analytic in C \ [0,i?]. 

(2) For large z, the resolvent has the asymptotic behaviour 

W{z) = ^ + ^ + 0{z~^) , |z|-)>oo, (2.7) 

where the leading coefficient is determined by the normalization condition 
(1^ . and E is the first moment of the equilibrium measure: 

rR 

MPo(Ai)dAi. (2.8) 

Jo 

(3) The particle density can be found as the discontinuity of the resolvent 
across its cuts, /q U /+, 

Po(M) = -^[W^(M + iO)-IT(M-iO)], (2.9) 

In particular, it follows that 

W{p + \Q)-W{p-\0)=r 

I —27ri p G 

(4) The resolvent satisfies 

W{p + iQ) + W{p + iQ) = V'{p), pGlo. (2.10) 

The above properties, together with some ‘ansatz’ on the structure of the partition 
of the interval [0,i?] into voids, bands, and saturation regions, are sufficient to 
completely determine the resolvent associated to the considered discrete Coulomb 
gas. 

Before discussing the issues related to the discrete nature of the Coulomb gas, 
and to the related phenomenon of emergence of saturated regions, it is useful to 
recall the procedure for the determination of the resolvent in a case where the 
discreteness and hard walls play no role, such as that of a single band [a,b], with 
two voids, [0,a] and [b,R]. The equation for the resolvent is a singular integral 
equation of the form (12.1011 . where Iq = [a,b]. The end-points of the band are 
determined by the condition that the solution of (12.101) should match at leading 
order the large z behaviour W{z) ^ 1/z, dictated by (12.711 . 
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However, if the density found by solving (12.101) does not satisfy the condition 
(1^ . or the end-points a and b are not internal to the interval [0, i?], then one 
has to introduce some suitable ansatz about the structure of the partition of the 
interval [0, i?] into voids, bands, and saturated regions, with consisting of at 
least one interval. 

The existence of at least one saturated region (that is, /+ ^ 0) implies the 
following form of the resolvent: 

W{z)= [ - dfi + H{z), H{z):= [ d^. 

Jl+ i/o ^ M 

Here the function H{z) can be found from the equation 

+ + = fielo, ( 2 . 11 ) 

where 

U{z) = -2 [ -^dn + V'{z). 

Ji+ z — 

Assuming that Iq consists of just a single interval (the so-called one-band ansatz, 
m, with the end-points a and b satisfying 0 < a < b < R, the solution of (j2.11l) 
such that H[z) = 0{\/z) as z —>■ oo, is 


H{z) 


^J{z-a){z-b) 

2-k 


U{p) 


{z - 


d/r. 


( 2 . 12 ) 


The end-points a and b can be found by solving the equations on the coefficients of 
order and z~^ terms of the large z expansion of (12.121) . to match the the large z 
behaviour (1^ at leading order. Finally, the explicit form of the particle density is 
derived by using (1^ . and its consistency against condition dUl) may be verified. 

In formulating some ansatz on the structure of the partition of the interval 
[0,i?] into voids, bands, and saturation regions, we take into account the form of 
the potential (12.41) . together with the presence of two hard walls (that is infinite 
potential barriers) at /r = 0, i?. It is clear that the band-gap structure may be 
affected by varying the parameters R G [l,oo) and Q G [0,oo). The whole solu¬ 
tion is obtained by investigating how the band-gap structure is modified as these 
parameters are varied. 

Concerning the shape of the potential (12.41) . it is clear that F(0) = 0 and, 
as /i — >■ oo, the function F(/i) is asymptotically linear, F(/i) ^ —(loga)/r, with a 
positive slope (recall that 0 < a < 1). The derivative of the potential reads 


V'ifJ.) = - log a -b log ^ - log(^ -b Q). 


Thus the function V{fj,) has a single minimum in the interval [0, i?], at /r = /tq, 
where 


A nonvanishing density (with bands and possibly saturated regions) is thus expected 
in the vicinity of fig- 

Concerning the presence of hard walls, it is known that these affect the band- 
gap structure in that they admit only gaps in their vicinity m- Thus, in our case, 
the intervals adjacent to the points /i = 0 and ^ = R can only be voids or saturated 
regions. It is clear, basing on the physical picture of the Coulomb gas, that the 
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presence of a void or a saturated region is related to the distance of the minimum 
Ho from the hard walls at /r = 0, i?. 

A detailed analysis provided below shows that four possible band-gap struc¬ 
tures, or scenarios may occur. We refer to them as Regime lA, Regime IB, Regime 
IIA, and Regime IIB for definiteness. In all these scenarios, there is just one band, 
that we denote [a, &], between two gaps, [0, a] and [b,R]. The four scenarios differ 
one another in the nature (void or saturated region) of the two gaps. Namely, in 
Regime lA the first gap, [0, a], is a saturated region, while the second gap, [b, i?], is 
a void; in Regime IB both gaps are voids; in Regime IIA both gaps are saturated 
regions; finally, in Regime IIB the first gap is a saturated region while the second is 
a void. The four regimes correspond to four different regions in the domain Q > 0, 
i? > 1. 


2.3. Regime I. Regime lA and Regime IB are closely related. Actually, they 
can be viewed as two different scenarios for one same regime. Regime I, character¬ 
ized by a void on the rigth gap [6, R]. The equilibrium measures for Regime lA and 
Regime IB were first derived in [13] . 

Proposition 2.1. The resolvent of the equilibrium measure for the continuum 
limit of the discrete Coulomb gas CH) in Regime lA is 


„rt X 1^1 V^a(z - b) + ^b{z - a) 

W{z) = - log Va - log A- V - 

yj{b- a)z 

_ + Q){z -b) + yCb -K Q){z - a) 

® ^J(b- a)z 

where the parameters a and b, the two band end-points, are 

_ (l — a(l -I- Q)) L _ \/^(l -I- Q)) 

1 — a ’ 1 — a 


(2.13) 


(2.14) 


In Regime lA, there is one band [a, 6], with a saturated region [0, a], and a void 
[6, R]. The saturated region in the interval [0, o] implies that the resolvent has the 
form 

W{z)=\og^—+H{z), H{z):= [ dfi- 

z-a Ja z- H 

The function Wif,^{z) should solve equation (12.1111 with U (z) = V'{z) and /q = [a, 6], 
where the function H{z) should solve (j2.11L with function U{h) chosen as 

U{h) = -loga-hlog-y—^. 

h{h + Q) 

Evaluating the integral in (12.121) with this function and coming back to function 
W{z), we obtain expression (12.1311 . 

The density in the interval [a, b\ is given by the expression 


Po{h) = ~ arctan 

TT 


a{b — p) 
5(/r — a) 


1 

H— arctan 

TT 


(a-l-Q)(b-H) 
(b + Q)(H-a)’ 


while, by construction, po(h) = 1, for p G [0,a], and Pq{p) = 0, for p G [b,R]. 
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Requiring W(z) ^ 1/z, as z ^ oo, from (j2.13l) we obtain the following equations 
for the end-points a and b: 


\/b -f Q — i/a -f Q !— -^(6 -l- Q){a + Q) + y/ab — Q 

-Tf- -p -= ^-= 1- (2-15) 

vo + va ^ 

These equations can be easily solved, with the result (|2.14|1 

It is clear from qualitative arguments that the current regime holds for relatively 
small values of Q, and relatively large values of i?, so that the minimum of the 
potential is relatively close to the left hard wall at origin, and relatively far from 
the right hard wall at R. 

Imposing a = 0, we get the value Q = Qc corresponding to the soft-edge/hard 
edge transition induced by the left hard wall at origin. We have 


Qc = 


I — a 


(2.16) 


It is clear that Regime lA corresponds to values of Q < Qc- 

Imposing b = R, we get the value R = Rc corresponding to the soft-edge/hard 
edge transition induced by the right hard wall located at R. We have 

(l -I- 


Rc = 


1 — a 


(2.17) 


Clearly, this is the minimum value of R for which the scenario of Regime lA is 
applicable. 


Proposition 2.2. The resolvent of the equilibrium measure for the eontinuum 
limit of the diserete Coulomb gas diu in Regime IB is 

yj a{z — b) + \/b{z — a) 


W{z) = — log y/a -I- log ■ 


(2.18) 


\/{a + Q)(z — b) -I- \/{b + Q){z — a) 
where the parameters a and b, the two band end-points, are given again by (|2.14p 


In Regime IB, there is one band [a, b] between two voids, [0, a] and [b, R]. The 
resolvent W{z) is simply given by the solution of (12.111) with U{z) = V'{z) and 
Iq = [a, b]. As a result, after evaluating the integral (12.1211 . we get expression 

(EH. 

The density in the band is given by 

. 1 / a{b — p) 1 I (a Q)(b — pi) . , 

/Oo(m =-arctanW--- - arctan W --, p. € [a,b], 

TT ^b(p,-a) TT ^ {b-\-Q)(^ — a) 

while, by construction, po{p) = 0, for p G [0, a] U \b, R]. 

Requiring W{z) ~ l/z, as z —>■ oo, from (|2.18l) we have the following equations 
for the end-points a and b\ 


\/b Q — \/a Q 

Vb — y/a 


= VQ! 


y^(a -I- Q){h + Q) - Vab- Q 
2 


Note that these equations differ from those in Regime lA, see (I2.15|) . just by the 
sign of yfa. Their formal solution, satisfying 0 < a < &, is given again by (12.141) , in 
which now Q > Qc, see (I2.16|) . 

We also note that the critical value of R, namely, its minimal value, for which 
the Regime IB is applicable, is given by (12.1711 . Thus the minimum value of R for 
which the scenario of Regime IB is applicable is the same as for Regime lA. 
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In Regime IB, all charges accumulate in the vicinity of the minimum and 
one should check that the obtained solution for (12.111) with U{z) = V'{z) and 
/q = [a,b] satisfies condition (12.61) . This is indeed the case, thus ensuring that no 
saturated region arises. 

Summarizing, a critical value Qc exists, that separates the two regimes, with 
Regime lA corresponding to Q < Qc and Regime IB to Q > Qc- The value Qc is 
determined by the condition a = 0, that is by the vanishing of the left gap, and 
the corresponding soft-edge/hard-edge transition. Similarly, a critical value i?c, 
determined by the condition R — b, exists, such that Regime lA and IB correspond 
to i? > i?c- The value Rc is determined by the condition b = R, that is by the 
vanishing of the right gap, and the corresponding soft-edge/hard-edge transition. 

2.4. Regime II. Regime IIA and Regime IIB are closely related. Actually, 
they can be viewed as two different scenarios for one same regime. Regime II, 
characterized by a saturated region on the rigth gap [6, R]. They both correspond 
to values R < Rc, with Rc given by (12.1711 . 

Proposition 2.3. The resolvent of the equilibrium measure for the continuum 
limit of the discrete Coulomb gas CH) in Regime IIA is 


W{z) = - log y/a + log 


z — R 


- log 


- log 


y/a{z — b) + y/b{z — a) 
y/iR- a){z- b) + y/{R- b){z- a) 
y/{a -\- Q){z — b) + y/{b -\- Q){z — a) 


(2.19) 


yJ{R- a){z - b) + yJ{R- b){z - a) 

where the parameters a and b, the two band end-points, are the solutions of the two 
equations: 

-y/R—b—y/R—a y/b-\-y/a 

^ ^ 

( 2 . 20 ) 


y/R — b + y/R — a y/b + Q — y/a Q 

+ V^(a + Q){b + Q)^ ^ _ fc) ^ p 

In Regime IIA there is a band in [a,b] between two saturated regions in [0,a] 
and \b, R] . This implies that the resolvent in this case has the form 


W{z) = log ■ 


log- 


-b 


+ H{z), H{z):= [ 

J a 


Poih-) 

Z — pi 


dpi. 


z — a z — R 

The function II{z) should solve equation (I2.11|) . where the function U{fi) is given 
by the expression 

TTi \ (M~o)^(^— 6)^ 

c/M = -logo + log 

Evaluating the integral (I2.12|) . and coming back to function W{z), we get (12.191) . 
The density is given by 


Po[p) = 1 H— arctan i 


l a{b- p) l_ 

b{pi — a) TT 


arctan ^ 


-arctan i 

TT 


' {a + Q){b- p) 
(b-\-Q){p-a) 

j{R-a){b- p) 

UR-b)ip-ay 


P e [a,b], 
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with po(m) = 1, for /X G [0, a] U [6, i?], by construction. 

Imposing bh(z) ~ 1/z, as z —>■ oo, from (j2.19p the system of equations (j2.20l) 
is obtained for the end-points a and b. 

Proposition 2.4. The resolvent of the equilibrium measure for the continuum 
limit of the discrete Coulomb gas (HH) in Regime IIB is 


W{z) = — log \/a — 2 log 


yj{b- a)[z - R) 


^{R-a){z - b) + ^/iR- b){z - a) 

\/a{z - b) + ^/b{z - a) 


log- 


( 2 . 21 ) 


y/(a+~Q)(z'--l^ + \/(b + Q)(z — a) 

where the parameters a and b, the two band end-points, are the solutions of the two 
equations: 

-y/R — b — \/R — a y/b—y/a 

( 2 . 22 ) 


y/R — b + y/R — a y/b + Q — yja -\- Q 
-^b+y/{a + Q){b + Q)-Q ^ ^ 

In Regime IIB there is a band in [a, &], with a saturated region in \b, i?]. There¬ 
fore the resolvent has the form 

Po(m) 


W{z) = \og^-^+H{z), H{z):= [ 

Z- R Ja 


Z — fJi 


■ d^. 


The function H{z) should solve equation (I2.11L where the function U{p,) reads 

CfM = -logon-log 

As a result, expression (12.211) is obtained for the resolvent. 

The density is given by 


PoifJ,) = — arctan n 

TT 


I a(b — p) 1 

b{p — a) TT 


-arctan n 


' {a Q){b - p.) 
ib + Q)ip-a) 


2 

-arctan 

TT 


{R — a) {b — p) 
{R — b){p — a)' 


p G [a,b], 


with po{p) = 0, for p G [0, a], and poip) = 1, for /x G [6, .R], by construction. 

The condition that the resolvent (12.2111 has the behaviour W(z) ^ l/z, as 
z —>■ 00 , leads to the system of equations (I2.22|) for the end-points a and b. 

Just as for the case of Regime I, Regime IIA and Regime IIB corresponds to 
values of Q in the ranges Q < Qc or Q > Qc, respectively. The value Qc is given 
by the condition that a = 0, that is by the vanishing of the left gap, and the 
corresponding soft-edge/hard-edge transition. In Regime II the critical value Qc 
has not the simple form (12.161) . but it is given instead by some more complicate 
expression involving the parameter i? in a non-trivial way. 

Finally, note that both in Regime IIA and Regime IIB, as the parameter R 
is decreased the band [a, 6] shrinks down, vanishing in the limit R —^ 1. More 
precisely, in Regime IIA, one has a —>■ & ^ 0, and in Regime IIB a,6—^-OasR—>-1. 
In both cases, the density tends to the saturated value 1 everywhere in the interval 
[ 0 , 1 ]. 
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2.5. The Arctic ellipse. Regime I and Regime II discussed above correspond 
to values of R in the interval R > Rc, and 1 < R < Rc, respectively, where 
i?c = Rc{Q) is given by (I2.17|) . and Q is allowed to take arbitrary values in the 
interval [0,oo). We discuss here how these conditions translate in terms of the 
coordinates x,y € [0,1]. 

Clearly, the transition between the two regimes corresponds to R = Rc{Q), 
that gives a relation between the two parameters R and Q, namely 

\/(l — oi)R = 1 + \/ 0(1 + Q). 

Recalling (|2.3L this relation translates into the following relation between the co¬ 
ordinates X and y: 

\/y = ^(1 - a)(l -x)- y/m. (2.23) 

This is readily recognized as the portion of the Arctic ellipes (11.91) , between the two 
contact points {x, y) = (1 — a, 0) and {x, y) = (0,1 — a). The additional condition 
x> y (corresponding to the condition i? > 0) selects the arc of ellipse A, see (11.101) , 
separating the two domains Vi and Vn, see Fig. H) 

It can also be easily seen that the values of x and y, which corresponds to 
R > Rc, i.e., to Regime I, are subject to the conditions y < yc and y < x, where 
Uc = yc{x) follows from ()2.23|) . and is given by 

_ [ (\/(l - ~ ~ X € [{I - ^/a)/2,l-a] 

^ |0 X G [1 — a, 1]. 

This conditions implies that in Regime I we have {x,y) G T>i. Similarly, one can 
verify that Regime II corresponds to values of the coordinates x and y such that 
(x,y) e Vu. 

Summarizing, Regime I of the Coulomb gas corresponds in the six-vertex model 
in an L-shaped domain to situations where the cut-off corner is sufficiently small to 
lie totally outside the Arctic ellipse. Regime II corresponds in the six-vertex model 
to situations where the cut-off corner is sufficiently large to overlap with the Arctic 
ellipse. As explained in the next section, this transition manifests as a third-order 
phase transition in the free energy of the Coulomb gas, and consequently, of the 
six-vertex model in an L-shaped domain. Instead, the transition between Regime 
lA and Regime IB (or also between Regime IIA and Regime IIB), although being 
as well a soft-edge/hard-edge transition, does not induce any discontinuity in the 
Coulomb gas free energy. 


3. The free energy 


To evaluate function $(i?, Q), we exploit the specific dependence of the discrete 
Coulomb gas (HH) on the variable a. Recall the definition ()2.8I) of the first moment 
of the equilibrium measure. Clearly, we have the relation 


E = lim 

s—^oo 


1 dlogIr^s,q 
d log a 


(3.1) 


Given the resolvent W{z), the quantity E can be easily extracted as the 
coefficient of its large z expansion, see (lO) . Using (ED in ED, we have the 
following ordinary first-order differential equation for the function $(i?): 


d 

9 log a 


$(R, Q) = E. 
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This equation can be easily solved, 

$(i?,Q) = 


(3.2) 


[ E— + C{R,Q). 

J a 

where C(i?, Q) is some function independent of a. 

To fix the integration constant C{R,Q), we use the fact that the function 
ip{x, y) admits direct evaluation when the parameter a takes the limiting values 
a = 0 and a = \. These two cases play the role of initial conditions in determining 
the function (p{x, y) in Regime I and Regime II, respectively. 

3.1. The limiting cases. Let us first consider the case a = 0. To start 
with, let us come back to the Coulomb gas at finite values of r, s, q. From (Hm) it 
follows that the expression Ir,s,q is a polynomial in a. Its lowest order term is of 
degree s(s — l)/2; the corresponding coefficient can be found (modulo the number 
of permutations of rrij’s) by setting mj = j — I, that gives 

Ir,q,s ^ Q, 

This, together with m, implies 

F I — 1 


Hence, by (II.Sp . we have, in particular. 

This condition will allow us to fix the integration constant C{R, Q) in Regime I. 

Let us now turn to the case a = 1. Since the function F’r.s.g has a zero at a = 1 
of degree s(s + q), the function ip{x,y) = ip{x,y]a) has a logarithmic singularity 
as a —>■ 1. Besides this, the function (p{x,y) has a regular a —>■ 1 part, which is 
relevant for our analysis. 

To compute it, we focus on the Hankel determinant appearing in F’r.s.g, see 
representation (HI). In the case a = 1, the determinant can be evaluated in closed 
form. 


Proposition 3.1. For arbitrary non-negative integers r,s,q, s <r, the follow¬ 
ing holds 


det 


E 

.m—O 


m- 


j+k -2 fm + q 


n 

j=0 


UKj + qY-^U + Q + rV- 

q\{r -j - l)!(2j + qy.{2j + g + 1)! ’ 


Indeed, the determinant in the left hand side can be recognized as the Gram 
determinant of Hahn polynomials Qj(rn;q,0^r — 1). The statement follows from 
the standard technique of orthogonal polynomials. For an account of the properties 
of Hahn polynomials, see, e.g., the monograph [31], Sect. 9.5. 

Therefore, as a —>■ 1, 


F, 


s-l 

n 


j'-U + <i)Kj + q + r)i 


1^1 + g)!(2j + g + 1)! 


• (1 


(3.4) 


Let s q = [xN], s = [yN], with x,y G [0,1] fixed, N = r s q. We set 

^s-l 


v^(x,j/)=- im ^logi n 


jKj + qY-U + q + rV- 


N—^oo 




g!(r -j - l)!(2j + g)!(2j + g + 1)! 
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(3.5) 


From (I3.4|) . recalling (11.81) . it follows that 

lim [(p{x,y) +xy\og{l - a)] = ^{x,y). 

oc—^l 

Using the Stirling formula, we have 

y) = log(a; + y) - ^ log(l - x - y) 


- — \ogx + 


(1 ^ n ^ y'^ ^ , (1 “ vf 


■log(l-a;) - y logy+ ^ 


log(l- 2 /). (3.6) 


Recalling (12.2L relation (13.51) implies 

= -{R + Q + lf'il^{x,y). 
Recall that x and y are related to R and Q by 

1 + Q 


1 


y = 


R + Q + l i? + Q + l 

see (12.31) . Finally, dne to (13.61) . we obtain the condition 

<3)L., =-^ logR + iog(R _ 1) _ (£±21! i„g(ij +(3) 

+ <1:!^ iog(i + Q) - <FL21! iog(2 + <3) 

+ + ( 3 . 7 ) 

This condition make it possible to fix the integration constant C{R, Q) in Regime 

II. 

After these preliminary remarks, we turn to the evaluation of the free energy 
of the discrete Coulomb gas in both regimes. 

3.2. Regime I. In order to extract the first moment of the equilibrium mea¬ 
sure in Regime IA and Regime IB, we evaluate the term 0(l/z^) in the z —>■ oo 
expansion of the resolvent, see (12.131) and (12.181) . respectively. In both regimes we 
obtain 

E = —-^■\/(a + Q){b + Q). 

Substituting here the expressions for the end-points a and &, see (I2.14L we have 

E = (Q + 1)--h -. 

1 — a 1 

Integrating with respect to a, and recalling dSSl), we get the free energy of the 
discrete Coulomb gas for Regime I, 

$(R, Q) = -(Q + 1) log(l - a) -I- log \fa + C{R, Q), 


where the integration constant C{R^ Q) is still to be determined. The last relation 
implies 

ipix,y) = C{R,Q), 

see p.2l) . Finally, the initial condition (13.3|) determines the value of the integration 
constant, 

C{R,Q) = 0. 
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As a result, we obtain 


^p{x,y) = 0, [x,y)&Vi, 


as stated in Theoreni ll.41 

3.3. Regime II: The end-point equations. In the previous section we ob¬ 
tained that the end-points a and b of the band in Regime IIA, see (12.201) . and 
Regime IIB, see (12.221) . obey the equations 



(3.8) 


vy/ah Y^(a -|- Q}(b -I- Q) — Q 
2 


where 



1 for Regime IIA 
1 for Regime IIB. 


These are essentially equivalent to quartic equations; a direct treatment is given in 
appendix A. The most important concern about the end-points in the approach we 
follow here, is that they need to be represented in such a way that the integral in 
(13.21) might be evaluated explicitly. 

For this reason we consider here a solution of equations (13.81) in which the 
second equation is viewed as defining an algebraic curve in the variables a and 
&, while the first equation is regarded as fixing a point on this curve. Indeed, by 
squaring the second equation twice (since it involves the linear combination of three 
square roots with different arguments), a and b are easily seen to lie on an algebraic 
curve of degree four. It turns out that in suitable variables, this is simply a cubic 
curve, which is moreover of rational (rather than elliptic) type. In this construction 
the first equation in (13.81) becomes a quartic equation for the parameter along this 
curve. Furthermore, the solution of this quartic equation can be viewed as a change 
of variables, so that the integration (j3.2l) turns out to be feasible. 

Instead of dealing with the end-points a and b as unknowns, we introduce the 
quantities 


(v^± lyy/aY 
4 


A± 


-f Q i V4" Q) 


(3.9) 


B± = 


4 


(\/ R — a i R — b^ 


C± 


4 


Obviously, the end-points can be expressed in terms of these new quantities, for 
example 


ct — A_|_ -)- A_ — 2-^/A_|_ A_, b — -)- A_ -)- 2-^A_. 


Similar expressions are valid also in terms of R’s and C’s. Namely, the quantities 
dSH) satisfy the multiplicative self-consistency relations 


A+A_ = B+B_ = C+C-, 


(3.10) 


together with the additive ones 


A+ + A_ = B+ + B_ - Q = R - C+ - C-. 
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(3.II) 

























(3.12) 


The first equation in (13.81) now reads 



and the second equation can be written in either of two forms 


A±+B±+ 2C± = N±, 


(3.13) 


where 


N+:=R + Q + 1, N_:=R-1. 


Obviously, (I3.10|) . (I3.11|) . (I3.12L and (13.131) constitute a system of six equations for 
the six unknowns A±, B±, and C±. 

Before turning to solution of this system, it is useful to mention that starting 
from one form of equation (13.131) . say, using the ‘+’ sign, the other form, with the 
’ sign, can be directly found from (13.111) . Indeed, rewriting (13.111) as 


A^ + A— + 0-1- + C— — R, 


(3.14) 


U-f + B— + 0-1- + 0_ — R Q, 


and summing these equations, it is clear that the two forms in equation (13.131) are 
equivalent. 

Dividing the two equations in (13.131) by \/A±B±, respectively, and using the 
relation A+/B+ = B-/A-, which follows from the first equality in (13.101) . a com¬ 
parison of the resulting relations yields 


N+ - 2C+ N_ - 20_ 


^A^ ^/AIbZ 


Therefore, using (13.101) . we can readily express A’s in terms of B's and O’s, as 
follows: 



(3.15) 




We may now use these expressions in (13.141) to eliminate the B's. 


Using (13.151) . and writing O’s in terms of z’s, from (13.141) we have 




Solving for U’s, we obtain 



(3.16) 


Recalling that B^B_ = C+C-, see (I3.10L we have 



Simplifying the z+z- term standing in both sides, we obtain the cubic equation: 


Q'^z+z- = {N-Z+ -I- N+z- — N+N+){z+ — z-)'^. 


(3.17) 
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Note that it contains only the third and second order terms in z’s, thus describing 
a rational curve. This means that z’s can be expressed as a rational function of 
some parameter along the curve. 

To obtain this parametrization, we set z+ = tz- in (13.171) and consider z- as a 
function of the parameter t. Then we get 


z_ 


1 

N-t + N+ 


■ QH 

At-ir 


+ N+N_ , 


z+ = tz-. 


Note that we need just that portion of this curve for which the C’s are positive. 
It can be easily seen that this corresponds to values of the parameter t in the 
interval [N+/N-,oo). It is useful to introduce a new parameter rj S [0,1], which 
parametrizes this portion of the curve as follows: 


N+l + r] 
N- 1 — r] 


In terms of the new parameter, z’s read 


2 \{N+{l + r,)-N.{l-r^)f j 


(3.18) 


As a result, recalling that C± = {E± — z±)l2, using expressions for B's (I3.16L 
and those for A’s (13.151) . we obtain that the parametrization (I3.18|) implies nice 
factorized expressions: 


{1 — ■q){l + Rri)(l + Q + {R + Q)ri) 

0+ = iV+- T -, 

(2 + Q + ( 2 i? + Q)tA 

^ A + ri){l + Q + Rri)[l + {R +Q)ri) 

(2 + Q + ( 2 i? + Q)ti') 

D AT A + v)A + Q + + Q + {R + Q)v) 

X3+ — iV-|_-T- 

(2 + Q T (2i? + Q)tA 
^ A-v)A + Rv)A + iR + Q)v) 

(2 + Q + { 2 R + Q)it) 

^ (1 + 'n){i + R'n){^ + {R + Q)-q) 

(2 + Q + { 2 R + Q)'r]) 

^ (1 — ??)(! + Q + Rr]){i + Q + {R + Q)r]) 

(2 + Q + (2i? + Q)t]) 


These expressions solve equations (13.101) . (13.111) . and (13.131) . and moreover they 
correspond to the solution in which all A’s, B’s, and C’s are positive, provided that 

??e[o, !]• 

Note that till now equation (I3.12p has never been used yet. To obtain the 
solution of the initial end-points problem, we exploit this last equation. Using 
(13.191) . we have 


^{1 + TjA{1 + Q + Rrj)(^1 + {R + Q)rA 

°‘{1-Tj)^{l + Rr]){l + Q + {R + Q)r]) ~ 
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(3.20) 


































This is a quartic equation, and we need only that root taking values in the interval 
[0,1]. It can be verified that such a root rj = r]{R,Q-,a) always exists, provided 
that the parameters R, Q, and a correspond to Regime II. An explicit expression 
for this root is not crucial for what follows; it is discussed in Appendix A. 


3.4. Regime II: The free energy. In order to extract the first moment 
of the equilibrium measure in Regime IIA and Regime IIB, we evaluate the term 
0(l/z^) in the z —>■ oo expansion of the resolvent, see (12.1911 and (12.2111 , respectively. 
In both regimes we obtain 

E = —-^■\/(a + Q){b + Q) + —y^{R- a){R - h). 

Expressing E in terms of quantities (13.9p . we have 

^ = T + 

To obtain the function $(i?, Q), we write the integral in (13.21) as 

<i>(i?,Q) = J E^^dr^ + C{R,Q), 

where the functions E = E{r]) and a = a{rj) follow from (13.1911 and (13.201) . respec¬ 
tively. 

Evaluating the integral, we find that 

$(R,Q) = fi!(i?,Q) + C(R,Q), (3.21) 

where the function n(i?, Q) = Q] rf) is given as a linear combination of seven 
logarithms: 


0(i?, Q)= Q + -R - log(l -V)+{^-R + Q- ^ log(l -I- 77 ) 

-f -f log(l -l- R'n) + ^2 ~ R E Q — RQ^ log(I Q Rt]) 

+ (^ - -r) log {1 + {R + Q)v) 

-l- i? -I- Q -l- RQ -\- log (1 -l- Q -l- (i? -I- Q)it) 

- log {2 + Q+{2R + Q)r,). (3.22) 

To find the integration constant C{R, Q), we recall that in our preliminary discus¬ 
sion we obtained the condition (13.71) . To employ this condition, we note that from 
(13.201) it follows that 77 —>■ 0 as a —>■ I, hence the constant C{R, Q) can be inferred 
by equating (13.211) at 77 = 0 with the expression (13.71) . that gives 


C(R,Q) 


n2 

logR- 




(1 -I- QY 


log(I -k Q) 


[R -l- Q -l-1)^ 


- log(R -k Q) 
log(i? -|- Q -|- 1 ). 


(3.23) 


Expressions (13.2111 . (I3.22p . and (I3.23p provide a solution of the free energy problem 
of the discrete Coulomb gas in Regime II. 
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From ( 1 ^ and (IQ) . it follows that in Regime II the function (p{x,y) is given 
by the following expression; 


y) = x{x, y; y) - xo{x, y), {x, y) e Pn, 


(3.24) 


where 


X{x,y]y) = -xy\og{l - a) + y log a 


+ (^y+a: 2 /-y-a;-y+^J log(l - v) + -2xy+^j log(l + y) 

(I - a; + l) ylog [y + (1 - x)y] - + xy - y - a;^ log [a; + (1 - x)y] 

- (I + a; - l) y log [y + (1 - y)y] - (^y - a;y + a;^ log [a; + (1 - y)y], 


(a; + yy 


log [a; + y + (2 - a; - y)y], (3.25) 


and the function Xoix,y), independent of y, reads 


X" , (1 — x) 


Xo{x,y) = —^loga; - 


r, ( 1 - 2 /)^ 


2 log(l - a;) - y logy- ^ 


log(l - y) 


(1 — a; — y)^ , , 

+ ^-log(l -x-y). 


Taking into account that the equation for the parameter y in terms of variables 
X and y reads 

^ (1 + yf [a; + (1 - a;)y] [y + (1 - y)y] ^ ^ 


and also using that 

1 — a = 


(1 - yY [y + (1 - a;)y] [a; + (1 - y)y\ 
y[x + y+{2-x- y)y]^ 


(3.26) 


(1 + yY [a; + (1 - a;)y] [y + (1 - y)y] ’ 

we may express a in terms of 77 in f|3.25p . thus obtaining for x(a:, y] 77 ) the following 
expression: 


(f — ^ ^ 

xix, y; y) = -xy log y H- - -log(l - y) + - log(l + y) 

- (1 - x)y log [y + (1 - a;)y] - a;(l - y) log [a; + (1 - y)y] 

+ a;(l - x) log [a: + (1 - x)y] + y(l - y) log [y + (1 - y)y] 

+ log [a; + y + (2 - X - y)y]. (3.27) 

We conclude by stating a remarkable property of representation (I3.24p . 
Proposition 3.2. Let: 

Then the following holds: 

Xo{x,y) = xix,y;h). 
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(3.28) 
















Verification is fairly straightforward. Representation (13.241) together with p.27l) 
and (13.281) immediately implies formula (11.111) . which is thus proven. 

The proofs of Propositions 11.71 and 11.81 are based on specific properties of the 
function y;??), namely that its first and second derivative with respect to 77 
vanish at r] = h. Further details and explicit calculations are given in Appendix B. 


Appendix A. End-point equations in Regime II: A direct treatment 

We discuss here a rather direct approach to treat the end-point equations (12.201) . 
Using notations (EH), we note that these equations can be written using only A_|_, 
B+, and C+, or only A_, R_, and C-, 

C± 


= iVa) 


±1 


A±+B±+ 2C± = N±. 

Introduce parameter 7 by 

A+ 

From the first and second way of writing of the end-point equations, we have, 
respectively, 

7 A+ ^ 


(A.l) 


A+ = 


A_ = 


l-|-7-|-2y^’ "" (1-I-7)^0-I-2,/7 

From the definition of the quantities A± we have 


— A_^ A_ — 2 A_|_ A_, b — -t- A_ -t- 2 A+A_ , 


\/ a + Q — . — Y^^p4^, 's/b + Q — + \/ 'jA-. 




and 


The expressions for the end-points a and b which follow from these two ways of 
writing are compatible with each other, provided that 

Q -)- A_|_ -(- A_ =-1- yA—. 

7 

Explicitly, this is a an equation for the parameter 7 , 

y/aN- . A+ 


Q+( 1 - 7 )- 


= (1 - 7 ): 


(1-I-7)y^-I-2.^7 '"1 + 7-1-2y^' 

Dividing by the expression standing in the right-hand side, 


Q 1 + 7 + 2y/(Ty ^ y/aN- 1 + 7 + 2,/drf 


N+ 1 - 7 ' A+ (1 + j)y/a + 2y/^ 

and introducing a new parameter w such that 


= 1 , 


costa = 


(l + 7)Va + 2^ 


1 + 7 + 2y/crf 
we may rewrite the equation in the form 

A B 


cos (jO sincj 
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sincj = 


= 1 , 


y/l - a{l - 7 ) 

1 + 7 + 2^/07 ’ 


(A.2) 


(A.3) 































where 


y/aN_ \/l -aQ 

iV+ 


(A.4) 


It is interesting to note that the three parameters (a, i?, and Q) in fact enters the 
equation only through their two combinations, A and B. 

It is possible to find the roots of (|A.3|) . by noting that it can be written as a 
quartic equation 


x'^ — 2Bx^ — (1 — — B^)x'^ + 2Bx — B^ = 0, x = sinw. 


The roots of this equation are 


B 1 
"= 2 + 2 ' 


^2 + 25'( 1 - cos ? ) ± + 25'( 1 - cos 


6» + 27r 


±a/b 2 + 25 I-cos 


' — 27r 


(A.5) 


where 


0 G [0 , tt), 6 = arccos ( 1 — 


2A^B^\ 

) 




and the signs are to be chosen independently. The root which lies in the interval 
[0,1] corresponds to the first and second signs in (IA.5P being ‘+’ and respec¬ 
tively, as it can be checked by expanding the expression in (jA.Sp in power series, 
say, in B, and noting that the required root has the behaviour x ^ {1 — A)~^B, as 
5 ^ 0 . 

To find a relation between ry and parameter (p, we first express parameter 7, 
dAll), in terms of 77 using (I3.19p . 

(1 -|- Rvi) (1 -b (i? -b Q)it) 

(1 -b Q -b Rri){l -b Q -b (i? -b Q)'ll) 

Using (13.201) we also have 


= 


(I - r]){l + Rr]) 


(1 -b?7)(l - Rri) 


(I -b (i? -b Q)r])[l -b Q -b (i? -b Q)ll) 


(I -b 77)(1 -b Q -b Rr]) ’ 
and substituting in the expression for, say, cosw, see (IA.2I) . we obtain 

A + rj 


cos UJ = \foL 


1-77' 


Consequently, 


77 = 


cos w — -y/a 
cos w -b \fa 


(A.6) 


An expression involving sin a; can be given by taking into account (IA.3I) . Thus 
(lA.sp . with the indicated choice of the signs, provides an explicit expression for 77. 
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Appendix B. Phase transition at the Arctic ellipse 

Here we outline some calculations related to the third-order phase transition 
at the curve A separating the regions Vj and Vji, see Fig. S] 

We first mention a particular parametrization of the points of A which appears 
useful in these calculations. Recalling that A is a portion of the ellipse 

(1-x-y)^ (x - yf ^ ^ 

a 1 — a ’ 

one can introduce an angle parameter cj), such that 

1 — X — y = \/a sin x — y = \/l — a cos 4>- 
If we further set a = sin^ A, A G (0, tt), then we have 

x = cos^ , 2 /= sin^ , {x,y) G A, (j) G [X,7r/2]. (B.l) 

To illustrate the convenience of this parametrization, let us consider our quartic 
equation in the case of {x,y) G A. We take it in its w-parametrized form (IA.3I) . 
Recalling relations (lO) . for the coefficients A and B, see (IA.4|) . we have A = 
•v/a(l — X — y) and B = \/l — Oi{x — y), and so (IA.3I) in this case reads 

sin^ A sin 6 cos^ A cos 6 

- 1 -;- = 1 - 

COS oj sm uj 

The root we need is simply w = tt/2 — </> (note that </> G [A,7r/2] and A G (0,7r/2), 
and hence w G [0,7r/2), as required). Substituting the obtained expression for the 
root in (|A.6|) . we get 

I sin — sin A /(/)-|-AA /(/) — AA I xy 

^\{x,y)€A sin ())-I-sin A *"°\2 / ^\2 J y {1 — x){l — y) 

In the last expression one can recognize the quantity h, see (|1.12|1 . and hence its 
interpretation as the value of ry on the Arctic ellipse. 

Let us now turn to study the properties of the function ip{x, y) = ip{x, y; a) in 
the vicinity of A. In particular, we want to verify that 

(p{x,y)C{a-acf, a ^ . 

Since Lp{x, y) = x(x, y\ ? 7 ) —x(a;, y\ h) where the function x(x, y, rj) is given by (13.271) . 
and since the dependence on a enters only through rj, the problem reduces to 
studying properties of x(a:, y; rf) as a function rj. By a direct calculation it can be 
verified that 

9yX{x,y;T])\^^f^ = 0, k = l,2. (B.2) 

At the same time the third-order derivative does not vanish, that gives a non-zero 
value of C = C{x,y), which can computed by the formula 


1 1 

1 1 

1 oi—OLc 

/a 1 AT 2/5 

[ [dr, log ay ' J 


26 































where the function a = airf) is determined by (13.261) . By a straightforward calcu¬ 
lation we obtain 

I 2 2 1 — X ^ — y 

V ^^^\n=h \ — h 1 + h~^ y + {1 — x)h a; -b (1 — y)h 

_ 1-a^ _ 1 - y ^ ^ (l-a:)(l-y) 

a: -b (1 — x)h 2 / -b (1 — y)h 1 — x — y 

However, in the case of i9j^x(a;, y; ?7) the corresponding calculation appears to be 
much more involved; the parametrization (IB.II) turns out to be very convenient (in 
particular, when resorting to computer symbolic calculations), providing the result 
in a nice factorized form: 

3 I (sin(() — sin A)(sin())-b sin A)^ 

=- 8»nAco»-Asin»» -' 

Taking into account that djj log a | = — (sin (/)-bsin A)^ / sin A sin </) and Uc = sin^ A, 

we finally obtain 

_ 2sin(i;ii-A) sinjcj} + X) _ ^x{l - x)y{l - y) 
sin'^ 2A 2ap(l — Oc)^ ’ 

and we recall that Oc as a function of x and y is given by (11.151) . 

Similarly, the third-order phase transition can be regarded as a property of 
the function (p{x, y) = ip{x, y\ a) with respect to the variables x and y, at fixed a. 
Namely, the function ip{x, y) = (p(x, y; a) has a discontinuity in its third normal 
derivative at the curve A, with vanishing first and second normal derivative. Let t 
be the variable along the normal of the curve A at the point (a;, y), then 

y)\ix,v)eA = 0, fc = 1,2. (B.3) 

Since ip(x, y) = x(a:, y; rf) — x(x, y; h), in computing the derivatives with respect to 
t one has to take into account only those terms which involve derivatives of the 
functions x(a:, j/; ij) and x(x, y, h) with respect to their parameters, y and h, respec¬ 
tively. Hence the property (IB.31) is just a consequence of the analogous property 
for the function x{x, y; y), see (IB.21) . For the third derivative we have 

dfy{x,y)\(^^yyA = [idty\r,=hf - (dthf] d^xix,y;y)\^^f,- 

Since dty\ri=h in general is not equal to dth, the third derivative is non-vanishing; 
its explicit expression can be obtained along the same lines as for the quantity C 
above. 
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